A-prophage induction modeled as a cooperative failure mode of lytic repression 
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We analyze a system-level model for lytic repression of A-phage in E. coli using reliability theory, 
showing that the repressor circuit comprises 4 redundant components whose failure mode is prophage 
induction. Our model reflects the specific biochemical mechanisms involved in regulation, including 
long-range cooperative binding, and its detailed predictions for prophage induction in E. coli under 
ultra-violet radiation are in good agreement with experimental data. 
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Viruses were one of the first biological systems to at- 
tract the attention of physicists Despite their appar- 
ent simplicity, it has become increasingly clear 0,0 that 
they provide a window into the full systems complexity of 
the cellQ, as well as representing a major evolutionary [B| 
and ecological force@ affecting all three domains of life. 
The viruses that infect bacteria are known as bacterio- 
phages (or phages). When some phages infect a bacterial 
cell, there can be two possible outcomes or pathways Q- 
In the lytic pathway, the phage hijacks the cell's machin- 
ery to replicate itself many times, and to release the repli- 
cates by breaking open or lysing the cell. In the lysogenic 
pathway, the phage integrates its genome into that of the 
host microbe, becoming a prophage, but otherwise does 
not damage the cell. This lysogenic state is very stableQ; 
however, an insult to the cell through, for example, star- 
vation or exposure to ultra-violet (UV) radiation[8( can 
trigger a process known as prophage induction [9(: the 
prophage is excised from the cell's genome, and viral 
replication occurs leading to cell lysis. The most well- 
studied lysogenic system is the bacteriophage A, or A- 
phage, which infects Escherichia coli. Understanding the 
lysis- lysogeny system in detail is important, because this 
system is one_of the simplest examples of a gene regula- 
tory network[l^| — a pervasive and fundamental form of 
biological organization and function, whose principles are 
still being elucidated. Although there has been consid- 
erable interest recently in the role of stochasticity ll| in 
the switching behavior betwe en ly tic and lysogenic states 
as part of the phage life-cycle 12|, |l3|, [l4|, |l5|, [l6| , here we 
focus on UV prophage induction, where a different mech- 
anism is involved. 

UV prophage induction experiments exhibit threshold 
behavior 17[, in which the fraction of induced lysogens, 
(i.e., prophage containing cells) rapidly increases as a 
function of UV dose. Under typical laboratory growth 
conditions, the fraction of induced lysogens versus the 
UV dosage obeys a power law with a power very close to 
4 [l?]]- Power law behaviors of this type can arise in sev- 
eral ways: (i) as an event caused by 4 independent hits on 



a "target" (target theory[18[); or (ii) a chemical equilib- 
rium reaction involving a substrate bound to 4 chemical 
species, and quantifie d by the empirical Hill equation for 
chemical kinetics[l9|, [2(|. Target theory and chemical 
kinetics could provide a way of understanding how UV 
dose curves can yield an exponent of 4, but have little 
connection to the biochemical regulatory mechanisms of 
A-phage lytic repression^, |21| . An alternative perspec- 
tive is to view UV induction in the framework of the 
standard stochastic model of lytic repression lj| 14 1 with 
adjusted rate constants. While these models are informed 
by the biochemistry, the mapping to UV prophage induc- 
tion remains dependent on unknown parameters. 

In this Letter, we show how the emerging understand- 
ing of the role of DNA loops and long-range cooperative 
binding 2^, in the biochemical picture of lytic repres- 
sion can account quantitatively for the phenomenology 
of prophage induction. Our approach is to abstract the 
biochemistry into a systems-level description, in which 
the lytic repressor circuit is represented as a device com- 
prised of a number of redundant elements and one fail- 
ure mode, lysis. This allows us to draw connections be- 
tween the biochemical regulatory mechanism and reliabil- 
ity theory[24|, |2^, [26| and also predict the characteristic 
power law for UV prophage induction. 

Biochemistry of lytic repression:- Fig. [T] illustrates a 
widely accepted model of the lytic repressor switch in A- 
phage (for a review see 0). The lytic repressor molecule 
CI dimerises and binds to specific DNA sites in the Or, 
and Or control boxes, OlX and Ol2, blocking expres- 
sion of genes under the control of the Pl promoter. In 
the Or control box, OrI and Or2 regulate Pr. Since 
only free OrI and OrI sites allow the expression of Pr 
controlled genes, CI binding at cither the OrI or Or2 
suffice for repression of Pr. Fig. [5] sketches the relation- 
ship between these 2 Or sites and Pr expression. The 
same applies to the role of OiX and 0^2 in suppressing 
Pl- Also, as drawn in Fig. [21 expression of genes under 
the control of both Pl and Pr promoters lead to lytic 
development of the prophage [30]. Derepression of all of 
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FIG. 1: Schematic diagram of the A-phage lytic repression 
system. Attachment of RNA polymerase to promoter regions, 
indicated here by arrows, lead to gene expression. Pr and Pl 
lead to expression of genes in the lytic pathway. However, CI 
dimer binding at OrI or Or2 blocks transcription of Pr while 
CI binding at OlI and Ol2 block transcription of Pl- Or3 
likewise regulates the transcription of genes by the promoter 
region Prm while CI bound to Or2 promotes transcription of 
genes from Prm- Dimers of CI are capable of forming stable 
quadramers when attached to adjacent sites such as at OrI 
and R 2^, @ (modeled by [27[|). Furthermore, Or\-Or2 
and Ol1-Ol2 quadramers can form a stable octamer in a 
long range interaction typically spanning 2.4 kb, but up to 
3.8kb[li[2i. 
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FIG. 2: A-phage lytic repression regulatory circuit. CI bound 
to OrI or Or2 blocks expression of lytic genes under the 
control of Pr. CI bound to OlI or Ol2 blocks expression 
of lytic genes under the control of Pl- Since both sets of 
genes under the control of Pr and Pl are required, blocking 
expression of either effectively represses lysis. Thus, each of 
the 4 CI bound components blocks lysis when intact. These 
can be seen as redundant elements that perform the same 
task. 



dose x, is 1 — pi for each of the 4 redundant components 
(the 4 CI dimers bound to OrI, Or2, OtX, and L 2) in 
the lytic repression system (see Fig. [2]). In general, the 
probability of failure as a function of total UV dose, F, 
for a system of n redundant components is 



the 4 binding sites results in lysis while bound CI dimers 
at any site block the lytic pathway. In UV induction, 
RecA-mediated autocleavage of CI monomers deprives 
the binding sites of available CI dimers. RecA-mediated 
autocleavage can happen once RecA is activated as part 
of the host SOS response to DNA damage 31 1. 

Abstraction of the lytic repression circuit:- Prophage 
induction can be understood as the failure of the lytic 
repression circuit, which consists of 4 redundant compo- 
nents that each prevent lysis. Each has a failure rate /ii 
(z = 1 ... 4) per UV dose x and a corresponding survival 
probability ^ = cxp(— \iiX). Fig. [3 diagrams the relevant 
aspects of the A-phage lytic repressor regulatory system. 
Each component consists of a lytic repressor CI dimer 
bound to one of 4 specific DNA sites, i.e., cither OrI, 
Or2, OiX, or O l 2. Each site regulates the expression 
of genes essential to the lytic pathway by its influence 
on either the promoter Pr, by OrI and OrI, or P L , by 
OiX and O l 2. Thus, these 4 components have redundant 
functionality, i.e., repressing lysis. Since suppression of 
genes under the control of either promoter keeps lysis 
in check, only damage to the final component results in 
lysis. 

Reliability theory of the repressor circuit:- To under- 
stand the failure rate of the system, i.e., the fraction of 
cells lysed, note that the probability of failure, for UV 



F = Y[(1- Pi ). 



(1) 



We model the effects of radiation on the failure rate of 
the lytic repression system, by assuming that near the 
threshold x c , pi = exp[— /i^(x — x c ) + 0{{x — x c ) 2 )]. 

By taking measurements of the fraction of failed sys- 
tems and Eq. [TJ the number of redundant elements in 
the system can be deduced. The fraction of failed units 
is then given by 

n 

F{x) = Y[{1 - exp[-iax]) » (px) n (fix < 1) (2) 
i=i 

where /i = (IX"=i Mi) 1 ^™ 1S an effective failure rate. The 
UV prophage induction curve describes the fraction of 
cells lysed as a function of UV dose x and is predicted 
to follow Eq. [3 with n = 4. In other words, the fraction 
of cells lysed can be computed from the effective failure 
rate /z for the 4 CI dimer bindings at OrI, Or2, OjX, 
and L 2. The effective rate of failure pL varies between 
different experimental systems and depends on a number 
of parameters. For example, different E. coli hosts may 
exhibit varying levels of RecA activity 31 , 12] • Alterna- 



tively, mutant cl alleles may offer operator site binding 
affinities 33]. Below, we test Eq. [2] against data from ex- 
periments on radiative induced lysis, and extract pi. 
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FIG. 3: The UV prophage induction curve. F/ (a: 4 ) versus UV 
dosage x. The data is taken from different A strains. Wild- 
type prophages were integrated as monolysogens and dilyso- 
gens with one or two inserted A genomes, respectively. All 4 
induction curves scale approximately as a power of 4 for small 
UV doses (0.2-2 J/m 2 ). (Inset) Log-log plot of the same data. 
Curves are right- and left-shifted so that they overlay each 
other. The line represents Eq. [T]for n = 4. 



Experimental tests:- In order to measure the frac- 
tion of lysogens induced as a function of UV dose, we 
followed standard protocols . Briefly, exponentially 
growing lysogenic cells were harvested and resuspended 
in buffer. The cells were then irradiated by a germici- 
dal UV lamp in dim ambient light for a range of doses 
at ~ 1 J/m 2 /s. After irradiation, aliquots were diluted 
into growth medium, shaken for 2 h at 37C in the dark, 
treated with CHCI3 and titered for plaque forming units. 

In our abstraction of the lytic repression system, we 
noted four redundant components, as shown in Fig. [51 
Since this implies n — 4, we can test the applicability 
of Eq. [2] to our experimental results by plotting F/(x 4 ) 
against UV dose shown in Fig. [3] 

In Fig. [31 the effective rate of failure p manifests itself 
as a horizontal line. As shown in Fig. [31 experimental 
data agrees with Eq. [5] for a certain range of UV dosage. 
Disagreement between theory and experiment occurs at 
both low and high UV dosage. The breakdown in both 
these regimes is readily interpreted. For extremely low 
dosage (near zero), spontaneous lysis events not induced 
by UV irradiation become the strongest contributing fac- 
tor to the failure rate. These failures come from other 
events, such as spontaneous RecA activity @ and muta- 
tions to A-phage|7j. At high radiation doses (not plot- 
ted, see [17(]), the fraction of lysis does not saturate at 
one and instead begins falling with respect to UV dose. 
Here, damage to the lytic pathway likely results in the in- 
ability of a cell to lyse, either because key components of 
the lytic pathway or host metabolism have been crucially 



damaged. In other words, this drop is not the effect of 
a lower failure rate of lytic repression, but a reflection of 
the high failure rate of other cellular systems upon which 
lysis relies. In these cases, failure of lytic repression can- 
not be detected by cell lysis since the lytic response has 
been disabled. 

Discussion:- Our model postulates that induction, or 
a failure event, is dominated primarily by CI dimer dis- 
sociation. Though many factors, such as DNA damage 
and RecA activity [35| , contribute in principle to the fail- 
ure of lytic repression, near the threshold of prophage 
induction, only the most rapidly-varying parameter is 
important. The rapidly- varying Cl-dimer operator site 
bindings about the threshold point for lytic induction x c 
justifies the approximation made in Eq. [5] allowing us to 
observe the predicted n = 4 power law. 

In contrast, chemical kinetic models assume that since 
the timescale of CI dimer dissociation (~ 30 sec[33j) falls 
an order of magnitude below the time required for lysis 
(~ 30 min), the CI bindings can be regarded as b eing adi- 
abatically slaved to CI monomer concentration 13. [T3|. 
and thus not a determining factor in the switch. The 
prophage induction curve and its power law behavior 
then arise from the behavior of the CI monomer deple- 
tion. However, as can be seen from the data presented 
in Fig. [31 the power law dependence n = 4 is robust, and 
not sensitive to the different strains or experimental con- 
ditions. In our model, the power n reflects the number of 
redundant elements — not the properties of the individ- 
ual components — and so is robust. This can be directly 
tested by manipulating the number of redundant opera- 
tor sites and measuring the resultant power law depen- 
dence. Ref. [13] measured Pr expression as a function of 
CI concentration. As shown in Fig. [1] Pr expression is 
regulated by 2 bound operator sites, so we predict that 
Pr expression should be described by a kinetic curve with 
an n — 2 power law dependence, as has been previously 
notedly. The kinetic data of ref. are indeed con- 
sistent with the prediction we have made here, based on 
our system-level abstraction of the underlying biology. 
This finding supports our view that the threshold behav- 
ior of prophage induction is determined by CI, and not 
by other steps in the repressor circuit. To establish this 
more conclusively, it would be necessary to check that Pl 
expression also exhibits the predicted n = 2 power law 
behavior. 

Our model implies that the repressor sites are the key 
to stabilizing the lysogenic state while the presence of Cro 
does not play a role in the switching 37|, [38( except to en- 
force commitment to the developmental transition[39j]. 
Once the switch has been activated, it cannot be re- 
versed, due to the role of Cro. Our model also suggests a 
mechanism for abortive induction events that are some- 
times observed [9]. In our model these arise when un- 
blocked Pl transcribes the genes required for excision 
(see [4j) while Pr remains blocked. 
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TABLE I: List of predicted and observed changes in inducibil- 
ity tabulated according to mutation, either in A or the E. coli 
host. | indicates that the change results an increase in fj,, 
while I indicates a decrease. Predictions based on reliability 
theory match with currently available data. 



Strain 



Phenotype 



Theory Data Ref. 



A cl ind s -\ faster CI cleavage i ] [32J 

lexA51 recA441 increased RecA activity [ j [32] 

A cl ind543 stronger CI dimerization J, J, [32] 

A cl Y210^N disrupts CI dimer-dimer j j [40] 
cooperative binding 

OrT weaker Or 2- CI binding ] ] [20] a 

"Measured Pr expression 

As shown by Fig. [31 different conditions lead to dif- 
ferent values for \x. We can use reliability theory to an- 
ticipate the trends in variation of \i between 2 similar 
experiments. The rate of component failure depends on 
a number of variables including RecA activity, CI con- 
centration, binding strength of repressor sites, and sta- 
bility of CI to autocleavage. Here, damage to redundant 
components corresponds to the dissociation of CI dimers 
from O l 1/2 and O r 1/2 sites (see Fig. [2J. 

Table [T] lists theoretical predictions for the variation in 
the failure rate of the lytic repressor or inducability aris- 
ing from possible laboratory manipulations of the rate 
of failure fj,. As expected, increasing RecA or rate of CI 
autocleavage leads to increased failure rates while an in- 
creased number of CI dimers slows down the rate of fail- 
ure. Increasing operator site binding strengths through 
cooperative binding or operator site mutations lowers the 
rate of failure. In summary, weakening the CI operator 
site bindings results in an increased failure rate for the 
lytic repressor, while increasing the probability of those 
binding events will have the opposite effect. 

We thank David Reynolds and Carl Woese for helpful 
discussions. IG wishes to express his deep gratitude to 
Edward Cox, at whose lab the phage induction exper- 
iments began. IG is partially supported by NIH grant 
R01-GM082837-01A1. NC was partially supported by 
Department of Energy grant No. DOE-2005-05818, and 
the IGB Postdoctoral Fellows Program. This material is 
based upon work supported in part by the National Sci- 
ence Foundation under Grant No. 082265, PFC: Center 
for the Physics of Living Cells. 



[1] E. Ellis and M. Delbruck, J. Gen. Physiol. 22, 365 (1939). 

S. Casjens, Mol. Microbiol. 49, 277 (2003). 
[3] N. Goldenfeld and C. Woese, Nature 445, 369 (2007). 
M. Ptashne, A Genetic Switch: Phage Lambda Revisited 
(Cold Spring Harbor Laboratory Press, New York, 2004). 
H. Brussow, C. Canchaya, and W. Hardt, Microbiol. Mol. 
Biol. Rev. 68, 560 (2004). 



[2] 
[3] 
M 

[5] 



[6] M. Weinbauer and F. Rassoulzadegan, Environ. Micro- 
biol. 6, 1 (2004). 

[7] K. Baek, S. Svenningsen, H. Eisen, K. Sneppen, and 
S. Brown, J. Mol. Biol. 334, 363 (2003). 

[8] S. Luria and R. Latarjet, J. Bacteriol. 53, 149 (1947). 

[9] A. Oppenheim, O. Kobiler, J. Stavans, and D. Court, 
Annu. Rev. Genet. 39, 409 (2005). 
[10] J. Hasty, D. McMillen, F. Isaacs, J. Collins, et al., Nat. 

Rev. Genet. 2, 268 (2001). 
[11] A. Arkin, J. Ross, and H. McAdams, Genetics 149, 1633 
(1998). 

[12] G. Ackers, A. Johnson, and M. Shea, Proc. Natl. Acad. 

Sci. U.S.A. 79, 1129 (1982). 
[13] M. Shea and G. Ackers, J Mol Biol 181, 211 (1985). 
[14] E. Aurell, S. Brown, J. Johanson, and K. Sneppen, Phys. 

Rev. E 65, 51914 (2002). 
[15] A. Lipshtat, A. Loinger, N. Balaban, and O. Biham, 

Phys. Rev. Lett. 96, 188101 (2006). 
[16] W. Bialek and S. Setayeshgar, Phys. Rev. Lett. 100, 

258101 (2008). 
[17] H. Kneser, Virology 28, 701 (1966). 

[18] G. Stent and J. Dohm, Molecular Biology of Bacterial 

Viruses (WH Freeman, San Francisco, 1963). 
[19] A. Hill, J. Physiol. (London) 40, 4 (1910). 
[20] N. Rosenfeld, J. Young, U. Alon, P. Swain, and 

M. Elowitz, Science 307, 1962 (2005). 
[21] D. Court, A. Oppenheim, and S. Adhya, J. Bacteriol. 

189, 298 (2007). 
[22] I. Dodd, K. Shearwin, A. Perkins, T. Burr, 

A. Hochschild, and J. Egan, Gene. Dev. 18, 344 

(2004). 

[23] I. Dodd, K. Shearwin, and K. Sneppen, J. Mol. Biol. 369, 
1200 (2007). 

[24] I. Gertsbakh, Adv. Appl. Probab. 16, 147 (1984). 
[25] L. Gavrilova and N. Gavrilova, J. Theor. Biol. 213, 527 
(2001). 

[26] L. Gavrilov and N. Gavrilova, in Handbook of Models 

for Human Aging, edited by P. Conn (Elsevier Academic 

Press, Burlington, MA, 2006), pp. 45-68. 
[27] J. Vilar and L. Saiz, Phys. Rev. Lett. 96, 238103 (2006). 
[28] B. Revet, B. von Wilcken-Bergmann, H. Bessert, 

A. Barker, and B. Muller-Hill, Curr. Biol 9, 151 (1999). 
[29] I. Dodd, A. Perkins, D. Tsemitsidis, and J. Egan, Gene. 

Dev. 15, 3013 (2001). 
[30] I. Dodd, K. Shearwin, and J. Egan, Current Opinion in 

Genetics & Development 15, 145 (2005). 
[31] J. Little and D. Mount, Cell 29, 11 (1982). 
[32] M. Dutreix, A. Bailone, and R. Devoret, J. Bacteriol. 

161, 1080 (1985). 
[33] H. Nelson and R. Sauer, Cell 42, 549 (1985). 
[34] J. Little, D. Shepley, and D. Wert, The EMBO Journal 

18, 4299 (1999). 
[35] J. Little, Proc. Natl. Acad. Sci. U.S.A. 102, 5310 (2005). 
[36] L. Bintu, N. Buchler, H. Garcia, U. Gerland, T. Hwa, 

J. Kondev, T. Kuhlman, and R. Phillips, Current Opin- 
ion in Genetics & Development 15, 125 (2005). 
[37] S. Svenningsen, N. Costantino, D. Court, and S. Adhya, 

Proc. Natl. Acad. Sci. U.S.A. 102, 4465 (2005). 
[38] S. Atsumi and J. Little, Proc. Natl. Acad. Sci. U.S.A. 

103, 4558 (2006). 
[39] R. Schubert, I. Dodd, J. Egan, and K. Shearwin, Gene. 

Dev. 21, 2461 (2007). 
[40] A. Babic and J. Little, Proc. Natl. Acad. Sci. U.S.A. 104, 

17741 (2007). 



